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O ! Abstract 

On 

We suggest a systematic procedure for classifying partial differential equations in- 

T S^ . variant with respect to low dimensional Lie algebras. This procedure is a proper 

synthesis of the infinitesimal Lie's method, technique of equivalence transformations 

^G ' and theory of classification of abstract low dimensional Lie algebras. As an application, 

& ■ we consider the problem of classifying heat conductivity equations in one variable with 

nonlinear convection and source terms. We have derived a complete classification of 
nonlinear equations of this type admitting nontrivial symmetry. It is shown that there 
are three, seven, twenty eight and twelve inequivalent classes of partial differential 
equations of the considered type that are invariant under the one-, two-, three- and 
four-dimensional Lie algebras, correspondingly. Furthermore, we prove that any par- 
tial differential equation belonging to the class under study and admitting symmetry 
group of the dimension higher than four is locally equivalent to a linear equation. This 
classification is compared to existing group classifications of nonlinear heat conductiv- 
ity equations and one of the conclusions is that all of them can be obtained within 
the framework of our approach. Furthermore, a number of new invariant equations 
are constructed which have rich symmetry properties and, therefore, may be used for 
mathematical modeling of, say, nonlinear heat transfer processes. 
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1 Introduction 

Traditionally group-theoretical, symmetry analysis of differential equations consists of two 
interrelated problems. The first one is finding the maximal Lie transformation (symmetry) 
group admitted by a given equation. The second problem is one of classifying differential 
equations that admit a prescribed symmetry group G. The principal tool for handling both 
problems is the classical infinitesimal routine developed by Sophus Lie (see, e.g., 0-0). 
It reduces the problem to finding the corresponding Lie symmetry algebra of infinitesimal 
operators whose coefficients are found as solutions of some over-determined system of linear 
partial differential equations (PDEs). 

Solving a classification problem for some group G provides us with an exhaustive descrip- 
tion of differential equations that are invariant with respect to this group and, consequently, 
could be analyzed by means of the powerful Lie group technique. And it is not just a matter 
of curiosity but the fundamental result that is used intensively in applications. An exper- 
imentalist, which believes that the nature is governed by symmetry laws, is provided with 
a criteria (symmetry selection principle) for choosing a proper nonlinear model describing a 
real process under investigation. Normally, a researcher has some freedom in choosing non- 
linearities of the model and it would be only natural to take those nonlinearities that provide 
the highest symmetry for the model. The classical example is the Lorentz-Poincare-Einstein 
relativity principle, which is to be respected by a physically meaningful model of relativistic 
field theory. From the point of view of the group theory the above principle is a requirement 
for a model under study to be invariant under the Poincare group (for more details, see, e.g., 
[0, [U ) . Consequently, finding all possible Poincare- invariant equations yields a complete ac- 
count of all possible ways to model processes of relativistic field theory by partial differential 
equations. 

In the overwhelming majority of papers devoted to solving classification problems a 
representation of symmetry group G (symmetry algebra g) is fixed. Given this condition, 
the problem is solved by a straightforward application of the Lie's algorithm. However, 
it becomes much more complicated if no specific representation of the symmetry algebra 
g is given. Then utilizing the Lie's algorithm directly one comes to the major difficulties 
arising from the necessity to find maximal symmetry algebra and solve classification problem 
simultaneously. A principal idea enabling to overcome the above difficulties was suggested by 
Sophus Lie. Indeed, his way for obtaining all ordinary differential equations in one variable 
admitting non-trivial symmetry algebras |5|, U teaches us what is to be done in the case in 
question. We should first construct all the possible inequivalent realizations of symmetry 
algebras within some class of Lie vector fields. If we will succeed in doing this, then symmetry 
algebras will be specified, so that we can apply directly the Lie's infinitesimal algorithm thus 



getting inequivalent classes of invariant equations. On this way, Sophus Lie has obtained his 
famous classification of realizations of all inequivalent complex Lie algebras on plane [||, [| . 
Recently, Lie's classification has been used by Olver and Heredero |7[] in order to obtain a 
classification of nonlinear wave equations in (1+1) dimensions that admit non-trivial spatial 
symmetries (i.e. symmetries not changing the temporal variable). What is more, Gonzalez- 
Lopez, Kamran and Olver || |9| have classified quasi-exactly solvable models on plane making 
use of their classification of real Lie algebras on plane || H . 

A systematic implementation of these ideas for PDEs has been suggested by Ovsjannikov 
[H]. His approach is based on the concept of equivalence group, which is the Lie transfor- 
mation group acting in the properly extended space of independent variables, functions and 
their derivatives and preserving the class of PDEs under study. It is possible to modify the 
Lie's algorithm in order to make it applicable for computing this group []]]. At the second 
step, the optimal system of subgroups of the equivalence group is constructed. The next 
step is utilizing the Lie's algorithm for obtaining specific PDEs belonging to the class under 
study and invariant with respect to the above mentioned subgroups. 

A further development of the Ovsjannikov's approach has been undertaken by Akha- 
tov, Gazizov and Ibragimov [|K], |11|] . They have obtained a number of classification results 
for nonlinear gas dynamics and diffusion equations. These ideas have been also utilized by 
Torrisi, Valenti and Tracina in order to perform preliminary group classification of some 
nonlinear diffusion and heat conductivity equations |12|, [13[]. Ibragimov and Torrisi have 
obtained a number of important results on group classification of nonlinear detonation equa- 



tions [1 i] and nonlinear hyperbolic type equations [15|. Note that there are number of papers 



(see, e.g., [L6| and the references therein) devoted to a direct computation of equivalence 
groups of some PDEs. Being somewhat more involved this approach has a merit of giving a 
possibility to find discrete equivalence groups or even non-local ones. 

The Ovsjannikov's approach works smoothly provided an equivalence group is finite- 
dimensional. However, if the class of PDEs under study contains arbitrary functions of 
several arguments, then it could well be that its equivalence group is infinite-parameter. The 
problem of subgroup classification of infinite-parameter Lie groups is completely open by now 
which makes problematic a direct application of the Ovsjannikov's approach. Consequently, 
there is an evident need for the latter to be modified to become applicable to the case of 
infinite-parameter equivalence groups. 

A possible way of modifying the Ovsjannikov's approach is suggested by the manner in 
which physicists construct nonlinear generalizations of the linear wave equations. They take 
a specific representation of the Poincare group realized on the solution set of the linear model 
and require that its nonlinear generalization should inherit this symmetry (for further details 
see, e.g., [HI). This approach makes the classification problem fairly easy to implement, since 



a representation of the symmetry algebra is fixed. A logical step forward is not to fix a priori 
a specific realization of the symmetry algebra but to fix the class of Lie vector fields within 
which this realization is searched for. It is namely this idea that enabled finding principally 
new nonlinear realizations of the Euclid [|J, Galilei [|, [L7|, |18| , extended Galilei [|17], |18| 



Schrodinger [|17], [18| , Poincare [|J and extended Poincare [19 , |2(| algebras. These results, in 
their turn, yield broad classes of Galilei- and Poincare-invariant nonlinear wave equations. 

What we suggest in the present paper is a proper combination of the above described 
approaches that enables a systematic treatment of a classification problem for the case 
of infinite-parameter equivalence group admitted by the class of PDEs under study. We 
perform group classification for the class of parabolic type equations describing nonlinear 
heat conductivity processes 

u t = u xx + F(t, x, u, u x ), (1.1) 

where u = u(t,x) is a smooth real-valued function, u t = du/dt, u x = du/dx and so on, 
F is a sufficiently smooth real-valued function. As shown below a direct application of the 
Ovsjannikov's approach is not possible since the equivalence group admitted by the above 
equation is infinite-parameter. By this very reason, a complete group classification has been 
obtained for particular cases of ( |1.1|) only []2T]]-|25 . 



The paper has the following structure. In the second Section we introduce the general 
method and necessary definitions and notions. The next section is devoted to computing and 
analyzing the equivalence group admitted by the class of PDEs ( |1.1| ). In Section 4 we carry 
out the preliminary group classification of ( |1 . 1|) , namely, we give a complete description of 
locally inequivalent PDEs of the form ( |1 . 1| ) that are invariant with respect to one-, two- 
and three-dimensional Lie algebras. In the fifth Section we present all inequivalent PDEs 
( |1 . 1|) admitting four- dimensional Lie algebras. Next, for each of thus obtained equations we 
compute the maximal Lie symmetry algebra thus obtaining the complete group classification 
of the corresponding models. In Section 6 we complete group classification of invariant heat 
conductivity equations with nonlinear source and show that there are no essentially nonlinear 
PDEs ( [Lip that admit symmetry algebras of the dimension higher than four. The seventh 
Section is devoted to an analysis of the connection of the results obtained in the paper to 
other classification results for ( |1 . 1|) known to us. It is shown that all of them can be derived 
from our classification of invariant PDEs ( |1 . 1[ ) . 

2 Description of the method. 

Our approach to group classification of PDEs is based on the following facts: 



PDE having a nontrivial symmetry admits some finite or infinite dimensional Lie al- 
gebra of infinitesimal operators whose type is completely determined by the structure 
constants. Furthermore, if the symmetry algebra is infinite dimensional, then it con- 
tains as a rule some finite dimensional Lie algebra (for example, the centerless Virasoro 
algebra contains the algebra sl(2, R).) 



Abstract Lie algebras of the dimension up to five have been already classified p6|, |27 



• Equivalence transformations preserving a class of PDEs under study do not change the 
structure constants of the Lie algebra admitted. 

Taking into account the above facts we formulate the following approach to group clas- 
sification of nonlinear heat conductivity equations ( |1 . 1|) : 



I. First of all we find the most general form of infinitesimal operators admitted by PDEs 
(|1 . 1|) . To this end we solve those determining equations that do not involve the function 
F . This yields a class X to which any symmetry of (|1 . 1| ) should belong. Next using 
infinitesimal or direct approach we construct the equivalence group Gs of the class of 
PDEs ( |1 . 1|) . Evidently, the group Gg sets an equivalence relation onl (two elements of 
Gg are called equivalent if they are transformed one into another with a transformation 
from Gg). We denote this relation as 8. 

II. At the second step, we find realizations of one-, two-, three-, four- and five-dimensional 
Lie algebras within the class I up to the equivalence relation £. To this end we use 
the classification of low dimensional abstract Lie algebras obtained by Mubarakzyanov 



26, 27 



III. Next, considering the obtained realizations of low dimensional Lie algebras as symme- 
try algebras of PDE ( |1 . 1| ) we classify all possible forms of functions F that provide 
invariance of the corresponding PDE with respect to this algebra. As a result, we 
get a complete classification of PDEs ( |1 . 1| ) admitting Lie symmetry algebras of the 



dimension up to five. 

IV. At the last step, we apply the Lie's infinitesimal algorithm for obtaining the maximal 
symmetry algebras admitted by those PDEs ( |1 . 1|) that are invariant with respect to 
four- and five-dimensional Lie algebras. This is being done straightforwardly, since the 
corresponding invariant PDEs (|1.1|) contains no arbitrary functions. 



Note that the above approach does not allow for a complete group classification of PDEs 



1.1|) , since there might exist realizations of higher symmetry algebras that does not contain 
four- or five-dimensional subalgebras. In fact, to get a full solution of classification problem 
one still has to be able to perform an exhaustive description of all inequivalent subalgebras of 
the Lie algebra of the infinite-parameter equivalence group Gg. However, in the case under 
consideration our approach enables solving the group classification problem for (|1 . 1|) in a full 



generality, since there are no essentially nonlinear PDEs of the form (|1.1|) whose symmetry 
algebra has a dimension higher than 4. 

It is also clear, how to modify the above approach in order to classify PDEs admitting 
some prescribed symmetry algebra (say, the Galilei algebra). At the second step, one has to 
fix the corresponding structure constants and find all inequivalent realizations of the Galilei 
algebra within the class X. Next, the maximal symmetry algebra is computed which yields 



the complete classification of Galilei-invariant PDEs of the form ( |I~T 



3 General analysis of symmetry properties of equation 

As a first step of group classification of PDE ( |1 . 1| ) , we find the most general form of the 
infinitesimal operator of the Lie transformation group admitted. Furthermore, we will con- 
struct the equivalence group of the class of PDEs ( |1.1| ). 

Following the general Lie's algorithm [jl], |] we are looking for an infinitesimal operator 
of the maximal symmetry group admitted by ( |1 . 1| ) in the form 

Q = rd t + td x + r ] d u , (3.1) 

where r = r(t,x,u), £ = £,(t,x,u), rj(t,x,u) are real-valued smooth functions defined in 
the space X ® U of independent t,x and dependent u = u(t,x) variables. The criterion for 
equation ( |1.1| ) to be invariant with respect to operator Q ( |3.1| ) reads as 



(p* _ <p** - rFt - £F X - V F U - tp x F Ux 
Here 



= 0. (3.2) 

(i.i) 



<p* = D t (r,)-u t D t (T)-u x D t {(i), 

<p x = D x {r,)-u t D x {T)-u x D x {£), (3.3) 

(p xx = D x ((p x ) - u tx D x (r) - u xx D x {£), 
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D t , D x are total differentiation operators defined in an appropriately prolonged space X®U: 

D t = d t + u t d u + u tt d Ut + u tx d Ux + . . . , (3.4) 

D x = d x + u x d u + u xx d Ux + u tx d Ut + .... 

Splitting ( |3.2| ) in a usual way and solving equations that do not involve F, we get the 
forms of the coefficients r, £ of the operator Q 

r = 2a(t), £ = a(t)x + b(t), 

where a(t), b(t) are arbitrary smooth functions and a(t) = %. Furthermore, the functions 
a(t), b(t), rj = f(t,x,u) and F(t,x,u,u x ) have to satisfy PDE 

ft - u x (ax + b) + (f u - 2a)F = f xx + 2u x f xu + u 2 J uu + 2aF t + (3.5) 

+ (ax + b)F x + fF u + f x F Ux + u x (f u - a)F Ux . 



Consequently, the maximal symmetry group admitted by equation (|1 . 1| ) is generated by 
an infinitesimal operator of the form 

Q = 2a(t)d t + (a(t)x + b(t))d x + f(t, x, u)d u , (3.6) 



functions a, b, f, F fulfilling the relation (|3.5|). 

Evidently, if we impose no restrictions on the choice of the function F, then the infinites- 
imal operator Q equals to zero and, consequently, the symmetry group of the nonlinear 
heat conductivity equation ( |1.1| ) reduces to a trivial group of the identity transformations. 
Non-trivial symmetry groups appears, if we specify in an appropriate way the source F. 

As we have mentioned in Introduction, there are different ways for constructing the 
equivalence group Gs for the class of PDEs ( |1 . 1|) . We use the direct method for finding the 
group G e - 

Let 

T — a(t,x,u), £ = f3(t,x,u), v = j(t, x,u) (3.7) 

be an invertible change of variables that transforms the class of PDEs fll.lj) into itself 

v T = Vtf + G(T,£,v,ve). (3.8) 

Computing the derivative u x yields 

v t ol x + v i fi x - 7 X . 



U.r 



lu - v T a u - v^(3 u ' 
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On the other hand, in view of arbitrariness of the function F it follows from ( |3.8| ) that the 
relation of the form 

u x = g(r,€,v,V() 

holds. Hence we conclude that in ( |3.7| ) a x = a u = 0, or a = a(t), a = ^ 7^ 0. 

Computing the derivatives u t , u xx with account of the relations a x = a u = <=> a = 
a(t), d^O we get 

u t = w T d(7 u -w 5 /5 n )" 1 + 6' 1 (r,^,w,^), 
u xx = v#{Pl(j u - v^f3 u )~ l + 2f3 x p u (v(:p x - j x ){j n - v^f3 u )~ 2 + 
+/3« W* - lx) 2 {lu - vs(3 u )~ 3 } + 9 2 (t, £, v, vt) 
with some function 6' 2 . Taking into consideration ( |3.8| ) yields the relation 

«(7« _ ^AJ 2 = /? 2 (7« - ^AJ 2 + 2f3 x f3 u {v^(3 x - 7 X )(7 U - ^/3 U ) + (3 2 u {v^(3 x - j x ) 2 . 

As a, 7,/? do not depend on u x , we can split the left-hand side of the above equation by v% 
thus getting the system of determining equations for the functions a, /3, 7 

(d - (3 x )-fl = y x P u (jxPu ~ Wxlu), 
-2(d - P 2 x h u p u = 2[3 2 x7u p u , 
a(3 2 u = 0. 

As d 7^ 0, it follows from the last equation that (3 U = 0. In view of this fact system in 
question reduces to a single equation 

(d - p 2 x H = 0. 

Since transformation of variables ( |3.7| ) is invertible, the relation 7« 7^ holds. Hence we 
get d = (3 2 . Consequently, d > 0, /3 = ±vda; + p(£). Summing up, we conclude that the 
equivalence group Gs of the class of PDEs (|1.1| ) reads as 



t = T(t), x = e\[t~(t)x + X{t), u = U(t,x,u), (3.9) 



where T(t) > 0, C/ u 7^ 0, T = g, e = ±1. 

Note that the infinitesimal method for finding the infinitesimal operator of the equivalence 
group yields the following class of operators (we skip the derivation of this formula): 

E = a{t)d t + -a(t)x + p(t) d x + r)(t, x, u)d u + [r) t - r) xx 

+ (Vu + Oi(t))F - u x (-a(t)x + p{t)) - 2u x r] xu - u 2 x r] uu }d F , (3.10) 



where a, p,r] — i](t, x, u) are arbitrary smooth functions. 

It is not difficult to become convinced of the fact that transformations ( |3.9|) can be 
obtained from the group transformations generated by operator ( |3.1(J| ) under condition that 



the latter is complemented by the discrete transformation x — > —x. Consequently, both the 
direct and infinitesimal approaches give the same equivalence group for the class of nonlinear 
heat conductivity equations (|1.1|). 

4 Preliminary group classification of equation ( ftTll ) 



In this section we classify equations of the form (|1.1| ) that admit invariance algebras of the 
dimension up to three. We start from describing equations admitting one- dimensional Lie 
algebras, then proceed to investigation of the ones invariant with respect to two-dimensional 
algebras. Using these results we describe PDEs ( |1.1| ) which admit three-dimensional Lie al- 
gebras. An intermediate problem which is being solved, while classifying invariant equations 
of the form ( |1 . 1| ) , is describing all possible realizations of one-, two- and three-dimensional 
Lie algebras by operators (|3.6| ) within the equivalence relation (|3.9| ). One more important 
remark is that PDEs that are equivalent to linear ones are excluded from further considera- 
tions. 

4.1 Nonlinear heat equations invariant under one-dimensional Lie 
algebras 

All inequivalent realizations of one-dimensional Lie algebras having the basis elements of the 
form Q3.6Q are given by the theorem below. 

Lemma 1 There are diffeomorphisms ( \3.9j ) that reduce operator ( \3.(\ ) to one of the following 
operators: 

Q = ±d t , (4.1) 

Q = d Xl (4.2) 

Q = d u . (4.3) 

Proof. Let an operator Q have the form (|3.6|). Making the transformation (|3.9|) we have 



Q^Q = 2aTd- t + 

+e(ax + b)\jT 



2a(X + ^xf(fy 



d x + [2aU t + {ax + b)U x + fU u ] fi . 



In a sequel, we have to differentiate between the cases / = and / ^ 0, that is why they 
are considered separately. 

Case 1. f = 0. Choosing U = U(u) in Q) yields 

Q = 2aT«9 f + [2a(X + -xf (T)~*) + e(dx + b)\ff]d s . 

If a = 0, then 6^0 (since otherwise the operator Q is equal to zero). So that choosing as 
T(t) in ( |3.9j ) a solution of the equation T = |6(t)| -2 we arrive at the operator 

Q = ±d s . 

Within the space reflection x — > —x we may choose Q in the form Q = d%. 

Given the inequality a 7^ 0, we put in ( |3.9| ) e — 1. Choosing as T(t),X(t) solutions of 
system of ordinary differential equations 

f -Wm = ^ 2a(t)X + 6(t)VT = 
2|a(t)| 

we arrive at the operator 

Q = ±dt- 



Case 2. f ^ 0. Provided a = b = 0, we can choose as U in (|3.9| ) a solution of PDE fU u = 1 
thus getting the operator 

Q = d u . 



If the inequality \a\ + \b\ 7^ holds, then choosing as U in ( p.9|) a solution of PDE 

2oC/ t + (ax + b)U x + fU u = 0, ^ ^ 

we come to the above considered case. 

It is straightforward to check that the operators ( |4.1|) - (|4.3|) cannot be transformed one 
into another with a change of variables ( |3.9|) . The lemma is proved. > 

Consequently, there are three inequivalent one-dimensional Lie algebras 
A\ = (ed t ) i Aj = (d x ), Al={d u ), e = ±l. 
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An easy calculation shows that the corresponding invariant equations from the class ( |1.1| ) 
have the form 



u t = u xx + F(x, u, u x ), (4.4) 

ut = Uxx + F(t,u,u x ), (4.5) 

u t = u xx + F(t,x,u x ). (4.6) 



A\ 

A\ 
A\ 

To proceed further, we need the transformations from equivalence group ( |3.9| ) preserving 
the forms of the basis operators of the above algebras. We give below the corresponding 
formulae 

A\ : t = £ + Ai, x = ex + A 2 , u = U(x 1 u) ) (4.7) 

A\ : i=t + X 1 , x = x + X(t), u = U(t,u), (4.8) 

A\ : t = T(t), x = e\frx + X(t), u = u + U(t,x), (4.9) 

{A l7 A 2 }cR, e = ±l. 

4.2 Nonlinear heat equations invariant under two-dimensional Lie 
algebras 

As is well-known, there are two different abstract two-dimensional Lie algebras, namely, 
the commutative Lie algebra A 2 .i = (Qi,Q2), [Qi,Q2] = and the solvable one A 2 .2 = 

(Qi,Q 2 >, [QiM=Q2. 



Theorem 1 The list of two-dimensional Lie algebras having the basis operators (3.6) and 



defined within the equivalence relation ( \3.Q ) is exhausted by the following algebras: 



A\ A = (d t ,d x ), A 2 21 = (d t ,d u ), 

A\ x = (d x , a(t)d x + d u ), A\ x = (d u , g(t, x)d u ),g ^ const, 

A\ x = (d x ,a(t)d x ), « = ^0; 

Al. 2 = {-tdt-\xd x d t ), Al 2 = (-2td t -xd x ,d x ), 

A 22 = \~uo u ,o u ), A 2 2 = {o x — uo u , o u ), 

A\. 2 = (ed t - ud u , d u ) , e = ±l. 
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Proof. Consider first the case of the commutative two-dimensional Lie algebra. Using Lemma 
1 we choose one of its basis operators (say, Q\) to be equal to one of those given in ( |4.1|) - 
Q4.3| ). For the sake of simplifying the form of the second basis operator Q 2 we make use of 
equivalence transformations (|4.7|) - (fL9|) . 

If Qi = ±<9i, then in view of the relation [Q\, Q2] = we obtain 

Q 2 = Xd x + f(x, u)d u , X = const. 



Provided the equation A = holds, taking as U in (|4.7| ) a solution of PDE fU u = 1 yields 
the realization A\ x . Given the inequality A 7^ we can choose as U in ( |4.7|) a solution of 
PDE XU X + fU u = 0, U u 7^ thus getting the realization A\ v 

Let us turn now to the case when Q\ = d x . Then the operator Q 2 takes necessarily the 
form 

Q 2 = Xd t + b(t)d x + f(t,u)d u , X = const. 



Provided A = 0, / 7^ 0, choosing as U in (|4.8|) a solution of PDE fU u = 1 we reduce the 
realization (Qi, Q 2 ) to become A\ v Next, if the inequality A ^0 holds, then taking as U, X 
in ( J4.8D solutions of system of PDEs 



AX + b = 0, XU t + fU u = 0, U u ^ 

we transform the operators Qi, Q 2 to the basis operators of the realization A\ l . The case 
A = / = gives rise to the realization A 21 . 

At last, consider the case when Qi = d u . Then 

Q 2 = 2a{t)d t + (ax + b)d x + f{t, x)d u . 

Utilizing the change of variables (|4.9| ), reduces the operators Qi, Q2 to the form 

Qi = du, 

Q 2 = 2afdt + [2a(e—=x + X)+}e\ff(ax + b)]d x 

2vT 
+ {2aU t + U x (dx + b) + f]d u . 

Given the conditions a = b = 0, we get the realization A\ 1 with / 7^ const. If a = 0, b 7^ 0, 



then choosing as T, U in (|4.9| ) solutions of system of PDEs 

/ T|6| = 1, bU x + f = 0, 
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we get the realization A\ A (a(t) = 0). 

Provided the inequality a^O holds, choosing as T, X, U in (|4.9|) solutions of system of 
PDEs 

2|o|T=l, 2aX + e^/fb = 0, 2aU t + U x (ax + b) + / = 0, 
transforms the operators Qi, Q2 to become 

Qi = du, Q2 = ±df 

thus yielding the realization A\ x . The fact that the obtained realizations of the two-dimen- 
sional commutative Lie algebra are inequivalent is established by a direct computation. 

Consider now the case of the solvable two-dimensional Lie algebra. Taking into account 
the results of Lemma 1 we analyze the three possible forms of the operator Q2 given in 

Let us first turn to the case Q2 = ±<9 t . In view of the automorphism of the algebra 
under study Q 2 — > —Q2 we may choose Q 2 = d t . Next, using the commutation relation 
[Qi,Q 2 ] =<?2 we get 

Q\ = (-£ + 2X)d t + (--x + S)d x + f(x, u)d u , A, 5 = const, 

where / is an arbitrary smooth function. 

Making use of the change of variables (|4.7|), where Ai = —2 A, A2 = —25 and U is a 



solution of PDE 

fU u + (5- l -x)V x = 0, U u ^ 0, 

we arrive at the realization A\ 2 . 

Consider now the case Q 2 = d x . Solving the commutation relation [Qi, Q 2 ] = Q2 yields 

Qx = (-2t + 2C 1 )d t + {-x + b{t))d x + fit, u)d u , d = const, 

where b, f are arbitrary smooth functions. 

Making the change of variables (|4.8|) with Ai = —C\ and X, U being solutions of system 
of PDEs 

2(C 1 -t)X + b(t)+X = 0, 
2{Ci-t)U t + fU u = 0, U u ^0 

transforms the operators Qi, Q2 to become 

Qt = -2idt - xd x , Q 2 = d x , 
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whence we get the realization A\ 2 . 

At last, consider the case Q 2 = d u . From the commutation relation [Qx, Q 2 ] = Q 2 we 
get the form of the operator Q x 

Q 1 = 2a(t)d t + (d(t)x + b(t))d x + (-u + f(t, x))d u , 



where a, b, f are arbitrary smooth functions. If a = b = 0, then choosing in (|4.9 ) U = —f 
we reduce the operators Q±, Q 2 to become 

Q\ = ~ud a , Q 2 = d^ 

thus getting the realization A| 2 . 

Provided a = 0, there exists a change of variables (|4.9| ) reducing the operators Qi, Q2 to 
the basis elements of the realization A\ 2 . The inequality a/0 gives rise to the realization 

^2.2- 

The fact that the realizations obtained are inequivalent is established by a direct verifi- 
cation. The theorem is proved. > 



Now we derive all inequivalent nonlinear heat conductivity equations (IT), that admit 
two-dimensional Lie algebras as symmetry algebras. 

For the realizations A\ x and A\ x the equations in question read as 

A\ x : Ut = u xx + F(u,u x ), (4.10) 

A 2 21 : u t = u xx + F(x,u x ), (4.11) 

corresp ondingly. 

Given the realization A\ x we may use the result of (|4.5| ) thus getting constraint (|3.5|) for 
the coefficient of the operator Q 2 in the form 

(xil x ± u . 

Hence it follows that 

F = —auu x + F(t, u x ) 

with an arbitrary smooth function F. 

So the most general PDE ( |1 . 1| ) invariant with respect to the Lie algebra A\ 1 reads 

A\ x : u t = u xx - auu x + F(t, u x ). (4.12) 

Treating the algebra A\ x in a similar way we represent constraint (|3.5|) as follows 

9t = 9xx + g x F Ux , g ^ const. 
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Given the relation g x = 0, the function g is constant, i.e., g = const. This means that 
PDE ( |1.1|) becomes linear. To avoid this we should impose the restriction g x ^ 0. Hence, 

F = (g t - 9xx)g^ 1 u x + F(t, x), g x ^ 0. 

Summing up, we conclude that the class of nonlinear PDEs of the form (|1.1|) invariant 
with respect to the algebra A\ 1 reads as 



/1 2.1 



Ut 



u , 



+ (9t- 9xx)g x 1 u x + F(t,x), g x ^o. 



(4.13) 



Turn now to the algebra A| x . Inserting the coefficients of the operator Q 2 into ( |3.5|) 
yields 

au x = 0, 

whence a = 0. This contradicts the assumption «^0. Consequently, there are no equations 
of the form (|1.1|) admitting A\ A as a symmetry algebra. 

Treating the algebras A\ 2 (i = 1, • • • , 5) in a similar way we get the following invariant 
equations: 



A 1 
^2.2 

A 2 
^2.2 

A 3 
^2.2 

^2.2 

^2.2 



u t = u xx + u 2 x F(u,xu x ); 
u t = u xx + t~ l F(u,tu 2 x ); 
u t = u xx + u x F(t,x); 
u t = u xx + u x F(t,e x u x ); 
u t = u xx + u x F(x, e et u x ), 



±1. 



(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 



Here F is an arbitrary smooth function. 

In what follows we will need equivalence transformations from the group Gs preserving 
the forms of the basis operators of all two-dimensional algebras considered above with an 
exception of the algebra A\ x . Omitting the derivation details we give the the subgroups 
of the group Gg that do no alter the forms of the basis operators listed in the assertion of 
Theorem 1. 



A 1 
A 2 
A 3 

A 1 
^2.2 



t = t + Ai, x = x + A 2 , u = U(u); 

t = t + Ai, x = ex + A 2 , u = u + U(x); 

t = t + Xi, x = x + X(t), u = u + U(t); 

t = T(t), x = e\/Tx + X( y t), u = u + U(t,x); 

t = t, x = ex, u = U(u); 



(4.19) 
(4.20) 
(4.21) 

(4.22) 
(4.23) 
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= t, x = x + X 1 Vt, u = U(u); (4.24) 

= T(t), x = e\ffx + X(t), u = u; (4.25) 

= t + Xi, x = x + X(t), u =u + e~ x U(t); (4.26) 

= t + Ai, x = ex + A 2 , w = u J r e~ t U{x)\ (4.27) 

Here {Ai, A 2 } C R, e = ±1. 

As the above transformations do not alter the form of the basis operators of the corre- 
sponding algebras, they can be used in order to simplify the form of the equations admitting 
the latter. An analysis shows that the only equation that can be simplified is PDE ( 4.13Q . 



A 2 

^2.2 


t 


A 3 

^2.2 


t 


A 4 

^2.2 


i 


^2.2 


i 



Indeed, the change of variables ( |4.22[ ), where T = t, X = and U is an arbitrary solution 
of PDE 

U t - U xx - (g t - g xx )g~ l U x + F(t, x) = 0, 

reduces Q4.13Q to the following equation (t — r, x — £, u — v ): 

v T = v# + (g T - 9tf)gz\, 

which is a particular case of fl4.16|) (up to notations). 

Thus equations ( |4.13| ), ( [4.1(j| ) are excluded from further considerations. 

4.3 Nonlinear heat equations invariant under three-dimensional 
Lie algebras 

We split the set of abstract three-dimensional Lie algebras into two classes. The first class 
contains those algebras which are direct sums of lower dimension ones. The remaining 
algebras are included into the second class. 

4.3.1 Equations ( |1.1|) invariant with respect to decomposable algebras 

The first class of Lie algebras contains two non-isomorphic algebras, namely, A3.1, A3.2. What 
is more, A 3A = (Qi,Q 2 ,Q 3 ), [Qi,Qj] = {i,j = 1,2,3), i.e., A 3A = A 1 ®A 1 ®A 1 = 3A t and 
A3.2 = (Qi, Q2, Qs), where [Q u Q 2 ] = Q 2 , [Qi, Qs] = [Q2, Qs] = 0, i.e., A 3 . 2 = A 2 . 2 © A x . 

Turn first to the case of the algebra A 31 . For describing inequivalent realizations of 
this algebra we use the results of the previous subsection on classification of inequivalent 
realizations of the algebra A 2 ,i, namely, of the realizations, A\ t , A\ 1 . 

Let A2.1 = A\ x . Then the relations Q\ = d t , Q 2 = d x hold, whence Q 3 = f{u)d u . Using 
transformation ( |4.19| ) we get the realization 

Qi=d u Q 2 = d x , Q 3 = d u . (4.28) 
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Consider next the case A 2 .\ = A\ v Then the relations Q\ = d t , Q 2 = d u , Q3 = \d x + 
f(x)d u , A G R hold. If A = 0, then we have the realization 

Qi = dt, Q2 = d u , Q 3 = f(x)d u , /V0. (4.29) 

If the inequality A 7^ holds, then using fl4.20|) with U being a solution of PDE XU x +f(x) = 
we come to conclusion that the operators Qi (i = 1, 2, 3) reduce to the form ( |4.28|) . 



Turn now to the case A 2 .\ = A\ A . In this case we have Q\ = d x , Q 2 = a(t)d x + d u , 
whence 

Q 3 = 2Ad 4 + b(t)d x + f(t)d u , 

where 2Ad = 0, A G R. If A ^ 0, then a = 0. Choosing as X, U in ( [4.2 1| ) solutions of 
system of PDEs 

2AX + 6 = 0, 2XU t + f = 

we reduce the operators Qi (i = 1,2,3) to the form (|4.28| ). Next, provided A = 0, the 
following realization is obtained 

Qi = d x , Q 2 = a(t)d x + d u , Q 3 = (3(t)d x + l(t)d u , (4.30) 

where a{t), j3{t),^{t) are arbitrary smooth functions such that the operators Qi ) Q 2 )Qz are 
linearly-independent . 

Thus, within the equivalence relations defined by ( |3.9| ), we have the three inequivalent 
realizations of the algebra A3.1, given by formulae ( |4.28|) -( |4.30|) Qi (i = 1,2,3). Now we 
proceed to constructing the corresponding invariant equations. 

Equation having as a symmetry algebra the Lie algebra ( |4.28| ) reads as 



u t = u xx + G(u x ), G Ux ^ const. 



The restriction for G guarantees that the above equation would not be of the form ( |4.16| ) . 



If the basis operators of the algebra A3.1 are given by ( [4.29 ), then F = F(x, u x ), and the 
invariance condition ( ^.5[ ) for the operator Q3 reads 

f" + fF Vm = 0, /VO. 

Hence it follows that 

F = -f"(fr 1 u x + G(x). 

As established above PDE ( |1.1|) with F = F, F being given by the above formula, is reduced 
to an equation of the form (|4.16| ) and therefore is not considered in a sequel. 
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Next, if the basis operators of the algebra A 3 ,± have the form ( |4.30| ), then F = —auu x + 



F(t, u x ) and what is more, the invariance condition ( |3.5|) for the operator Q 3 takes the form 

7 = - ja)u x , 

then 7 = Ci, (3 = 7a + C2, {Ci, C2} C R. In view of this fact, we have Q 3 = = C\(ad x + 
d u ) + Cidx = C1Q2 + C2Q1, which contradicts to the requirement of linear independence of 
the operators Qi {i = 1, 2, 3). 

Summing up, we conclude that there is only one realization of the algebra A 3 i, which is a 
symmetry algebra of PDE belonging to the class (|1 . 1| ) and cannot be reduced to an equation 
of the form ( f4.16| ) . Namely, we have 

^3.1 = (9t, d x , d u ), 

ut = Uxx +G(u x ), G Ux 7^ const. 

Let us turn now to analysis of realizations of the algebra A 32 = ^2.2 © Al- In order to 
describe these we use the realizations A\ 2 , ^22' ^22' ^2 2 °f the two-dimensional algebra 
A2.2 obtained in the previous subsection. 

Consider first the case when A2.2 = A\ 2 . Then Q\ = —td t — \xd x , Q 2 = d t , Q 3 = 



f(u)d u , f 7^ 0. It is not difficult to check that transformation ( |4 . 2 3| ) , where U is a solution 
of PDE fU u = 1, reduces this triplet of operators to the form 

Qi = -td t - -xd x , Q 2 = d u Q 3 = d u . (4.31) 

Next we turn to the case when A2.2 = A\ 2 . With this choice of A2.2 we get Q\ = 
—2td t — xd x , Q 2 = d x , Q 3 = XJ\t\d x + f(u)d u , A G R. If A = 0, then / 7^ and we arrive 
at the realization 

Q 1 = -2td t -xd x , Q 2 = d x , Q 3 = d u . (4.32) 

Provided / = 0, A 7^ 0, we have the realization 



Q x = -2td t -xd x , Q 2 = d x , Q 3 = ^\t\d x . (4.33) 

At last, if the inequality Xf 7^ 0, holds, then within the transformation ([4.24) we obtain the 



realization 

Q x = -2t8 t - xd x , Q 2 = d x , Q 3 = J\t\d x + d u . (4.34) 

The case A2.2 = ^22 gives rise to the realization 

Qi = d x - ud u , Q 2 = d u , Q 3 = Xd t + b(t)d x + e~ x f(t)d u , X e R. 



If A = b = 0, then the following realization is obtained 

Q x = d x -ud ui Q 2 = d u , Q 3 = e~ x f(t)d u , / ^ 0. (4.35) 

Next, given the conditions A = 0, b ^ 0, we can choose in fl4.26| ) U = b~ l f and reduce the 
initial operators to the form 

Q l = d x - u d u , Q 2 = d u , Q 3 = a(t)d x , «/0. (4.36) 

If the inequality Xb ^ holds, then we arrive at the realization 

Qi = dx - ud u , Q 2 = d u , Q 3 = d t . (4.37) 

Consider next the case when A 2 . 2 = ^ 22 - Then we have 

Qi = ed t - ud u: Q 2 = d u , Q 3 = C x d t + C 2 d x + e~ et f(x)d u , 

where {C*i, C 2 } C R, e = ±1. Hence we get within transformations ( |4.27|) and the choice of 
the basis the following three realizations: 

Q x = ed t - ud u , Q 2 = d u , Q 3 = e- et f(x)d u , f ^ 0, e = ±1, (4.38) 
Qi = ed t - ud u , Q 2 = d u , Q 3 = d x , e = ±1, (4.39) 

Q t = edt-ud u , Q 2 = 8 U , Q 3 = d t + Xd x , A > 0, e = ±1. (4.40) 

Evidently, the above obtained realizations of the algebra A 32 and the realization A\ ± are 
inequivalent. 

Now we choose from the set of so obtained realizations of three-dimensional Lie algebras 
those which are subalgebras of symmetry algebras of PDEs ( |1 . 1| ) not reducible to the form 

(ED- 



Equation invariant with respect to the algebra (|4 . 3 1|) reads as 

Ut = u xx + u x G(u), us = xu 2 x , G 7^ Xu^ 1 , A G R. 

Similarly, we get PDE of the form (|1 . 1| ) admitting the algebra A 3 . 2 having the basis 
operators fl4.32| ) 

ut = Uxx + t~ l G(uj), u = tu x , G 7^ X^/cJ, A G R. 

If we have realization ( f4.33| ), then F = t~ x F(u, tu%) in ( |1.1| ). That is why invariance 
condition (|3.5| ) for the operator Q 3 takes the form 

Ut = 0, e = ±1. 



2\/\t\ 
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Hence we conclude that there are no PDEs of the form ( |1 . 1| ) invariant with respect to the 
algebra under consideration. 

Provided we have realization (PI), the function F takes the form F = t~ 1 F(u, tu 2 x ) and 



invariance condition Q3.5|) for the operator Q 3 reads as 

2J\i 



where e = 1 under t > and e = — 1 under t < 0. Consequently, 

F = — -J \uj\u + G(uj), lu = tu x 

and the invariant PDE is given by the following formula: 

1 



Ut = u xx — ^-t 1 uJ\ui\+t l G(u), on = tu 2 x , G ^ \\J\uj\, A G R. 



Now we turn to the case when the operators Qi (i = 1, 2, 3) take one of the forms ( |4.35| ) 



(|4.37|) . If this is the case, then F = u x F(t,e x u x ) and invariance condition ( |3.5|) for the 



operator Q 3 is given by one of the corresponding formulae below 

/ = f(l-F-uF u ), uo = e x u xl 
—a = auFu, oj = e x u x , 
F t = 0. 

Integrating these PDEs yields the forms of the functions F in ( |1 . 1|) 

F = u x (ff- l -l) + e x G(t), 

F = — dtoT u x \w \uj\ + u x G(t), uj = e x u x , 

F = u x G(oj), oj = e x u x , G ^ Acj -1 , A e R. 

A further analysis shows that only the second and the third expressions for F from the 
above list give rise to essentially new PDEs of the form (|1 . 1| ) . 

At last, similar reasonings for triplets (|4.38|) - (|4.40|) give the following expressions for the 
functions F 

F = -(/ + /")(/'K + e- £f G(x), 

F = u x G(e €t u x ), G 7^ const, 

F = u x G(uj), uj = (u x ) x e e{xt - x \ A>0, G ^ const, e = ±1. 
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Again, only the second and the third expressions for F from the above list give rise to 
essentially new PDEs of the form fll.l|) . 

We summarize the results on classification of nonlinear heat conductivity equations ( |1.1[ ) 
invariant under the three-dimensional Lie algebras belonging to the first class in Table 1, 
where we use the following notations: 

(d t ,d x ,d u ); 
-td t - -xd x ,d t ,d u ); 



and what is more, e = ±1 



4 1 - / 

-^3.1 — \ 


4 1 - / 

^3.2 — \ 


4 2 / 

^3.2 — \ 


4 3 / 

^3.2 — \ 


4 4 / 


4 5 - / 


4 6 - / 

^3.2 — \ 


4 7 - / 

^3.2 — \ 


= ±1. 



-2td t -xd x ,d x ,d u ); 

-2td t - xd x , d x , \J\t\d x + d u ); 

d x -ud u ,d u ,a(t)d x ), «^0; 

d x -ud u ,d u ,d t ); 

ed t -ud u ,d u ,d x ); 

ed t - ud u , d u , d t + \d x ), A > 0), 



4.3.2 Equations ( |1.1|) invariant with respect to non-decomposable algebras 

Here we consider those three-dimensional real Lie algebras A 3 = (Qi, Q2, Q3) that cannot be 
decomposed into a direct sum of lower dimensional Lie algebras. The list of these algebras 
is exhausted by the two semi-simple Lie algebras 

A 3 . 3 : [Qi,Q 3 ] = -2Q 2 , [Qi,Q 2 ]=Qi, [Q 2 ,Q 3 ] = Q 3 ; 
A 3A : [Qi,g 2 ] = Q 3 , [Q2,Qs] = Qi, [Q 3 ,Qi]=Q 2 ; 

nilpotent Lie algebra 

^3.5 : [Q2, Q 3 ] = Qi, [Qi, Q2] = [Qi, Q 3 ] = 

and six solvable Lie algebras (non-zero commutation relations are given only) 



^3.6 
A3.7 
A3.8 
^3.9 
^3.10 
A3.11 



[QiM = Qi, [Q2,Q 3 ] 

[Qi,Qs] = Qi, [Q2,Qs] 

[QuQ 3 ] = Qu [Q2,Q 3 ] 

[QiM = Qu [Q2,Q 3 ] 



Q1 + Q2; 

Q2; 

-Q 2 ; 

qQ 2 (0 < \q\ < 1); 



[Qi,Q s ] = -Q2, [Q2,Qs} = Qi; 

[Qi, Qs] = qQi - Q2, [Q2, Qs] = Qi + <?Q 2 , g > 0. 
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Table 1. Equations ( |1.1| ) admitting the algebras A3.1, A 3 . 2 



Algebra 


Function F 


A 1 


G(u x ), G Ux ^\ AgR 


A 1 


u x G{ui), u = xu x , G ^ \uj~ 2 , A G R 


A 2 


r 1 ^^), UJ = tu 2 x , G^Xyfc, AGR 


^3.2 


— |t _1 Mi/cJ +t _1 G(cj), CJ = tU x 


A 3.2 


—aa~ l u x In cj + u x G(t), «^0,w = e x « x 


A 5 
^a.2 


Wa;(j(a;), oj = e x u x , G ^ Au; -1 , A G R 


^3.2 


u x G(u), u) = e et u x , G 7^ Xu^ 1 , A G R, e = ±1 


A 7 


u x G(u), uj = (u x ) x e< xt - x \ A > 0, G ^ const, e = ±1 



While constructing inequivalent realizations of the above algebras within the class of 
operators (|3.6|) , we use wherever possible the classification results obtained for the lower 
dimensional Lie algebras. 

Consider first the semi-simple algebras. Let A 3 = A 33 . Then Qx, Q 2 satisfy the commu- 
tation relation [Qx, Q 2 ] = Qx and form a basis of a two-dimensional Lie algebra isomorphic 
to A 2 .2- Indeed, choosing Q 1 = Q 2 Q 2 = —Q x we see that [Q x , Q' 2 ] = -[Q 2 , Qi] = Qi = Q' 2 - 
Thus we can use the results on classification of the algebra A 2 . 2 . According to the results of 
Subsection 4.2 studying realizations of the algebra A3.3 reduces to finding the form of the 
operator Q 3 for each pair of the operators Qx, Q 2 given below 



1) Qx = ed t , Q 2 = td t + -xd x ; 

2) Qx = d x , Q 2 = 2td t + xd x ; 



(4.41] 
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3) Q\ = d u , Q 2 = ud u - d x ; 

4) Qi = d u , Q 2 = ud u - ed t . 

Here e = ±1. 

One more remark is that the form of the operator Q 3 can be simplified with the use of 
transformations flOg) , (gg| ), (ggg ), ( ^2p . 

Let Qi, Q2 be given by the first formula from ( |4 . 4 1 ) . Then it follows from the commuta- 
tion relations 

[Qi,Q 3 ] = -2Q 2 , [Q 2 ,Q 3 ] = Q3 (4.42) 

that 

Q 3 = -et 2 d t - etxd x + x 2 f{u)d u . 

Given the condition f(u) 7^ 0, change of variables (|4.23| ) with e = 1 and U being a solution 
of PDE fU u = 1 reduces the operator Q3 to the form 

Qs = -ePdi - eixd s + xdy,. 

Consequently, we get the realization 



Qi = d t , Q 2 = td t + -xd x , Q 3 = -t 2 d t - txd x + ex 2 d u , s = 0, 1. 



(4.43) 



Let Qi, Q 2 be given by the second formula from ( |4.41|) . Checking commutation relations 
( |4. 42 ) yields that there is no operator Q 3 of the form ( |3.ti| ) which enables extending the 
algebra A 2 . 2 to the algebra A 33 . The same assertion holds for the remaining pairs of operators 
from QL41| ). 

Thus there exists a unique realization of the algebra A 33 that is given by ( f4.43| ). In this 
case, F = u 2 x F{u,u)), to = xu x and, consequently, invariance condition ( |3.5|) for the operator 
Q 3 takes the form 

e^Fu + 2uj 2 F u + AuF + 2] = -ecu. 

Provided E\ = 0, we get the equality uo = whence it follows that the only possible value of 
£1 is £1 — 1. With this condition, 

F = -- - u)- 1 + l)- 2 G(2u - u). 
4 

Hence we conclude that equation (|1.1|) invariant with respect to the algebra 



^3.3 = ( ed t, td t + ~xd x , -et 2 d t - etxd x + x 2 d u ) 
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reads as 

Ut = u xx + £-u x — x~ u x + x~ 2 G(2u — xu x ), e = ±1. 

On having used the equivalence transformation 

t — > t, x — > x, u — > —u, 

we may choose e = 1. 

Note that the algebra A\ 3 is isomorphic to the Lie algebra of pseudo-orthogonal group 
0(1,2). 

Turn now to the algebra A 3A . It does not contain a two-dimensional subalgebra and we 
use the classification results for one-dimensional algebras (Subsection 4.1). According to 
these results the operator Qi is reduced to one of the following inequivalent forms 

±d t , d x , d u . (4.44) 

Given the relation Qi = ±<9 t , we verify that there are no operators Q 2 iQz of the form 
T0[) satisfying together with Qi the commutation relations 

[Qi, Q 2 ] = Qs, [Q 2 , Qs] = Qi, [Qs, Qx] = Q 2 - 



Consequently, the class of operators ( PTB| ) does not contain operators Q 2 , Q3 that extend 
a realization of the one-dimensional algebra (Qi) to a realization of the algebra A 3A . The 
same assertion holds true for the remaining realizations of the operator Qi, Summing up 



we conclude that there is no PDE of the form ( |1 . 1| ) whose symmetry algebra contains a 
three-dimensional algebra isomorphic to A3.4. 

The algebra A3.5 contains the commuting subalgebra having the basis operators Qi,Q2- 
Since the latter is isomorphic to the Lie algebra A 2 .i, we can use the results of Subsection 
4.2. In view of these we conclude that there are three inequivalent realizations of the algebra 
A2.1 which might be invariance algebras of equations of the form ( |1 . 1|) , namely, 

A\. x = (d t ,d x ); 

A\ A = (d t ,d u ); (4.45) 

A\ A = (d x ,a{t)d x + d u ). 

Therefore, while considering the algebra A 35 we can suppose that Qi,Q 2 are given by 
one of the formulae ( f4.45|) . In order to simplify the form of the operator Q3 we use transfor- 



mations (fCTg), QODQ , QOT]) , respectively 
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Let the operators Qi,Q2 form a basis of the algebra A\ v If Q x = d x , Q2 = d t , then 
analyzing the commutation relations 

[QuQa] = a, [Q2,Q 3 ] = Qi (4.46) 

yields that the class of operators (|3.6|) does not contain an operator Q 3 which forms together 
with Q\,Q2 a basis of the algebra A 3 _ 5 . 

Next, provided Qi = d x , Q2 = d t , it follows from (|4.46|) that 

Qs = (t + X 2 )d x + f(u)d u . 
There is a transformation ( |4. 18| ) that reduce Q 3 to the form 

Q 3 = td x + ed u , e = 0, 1. (4.47) 

The most general PDE ( |1 . 1| ) , which is invariant with respect to the algebra A\ x reads 

u t = u xx + F(u,u x ). (4.48) 

That is why, condition for PDE ( [L.lfl to be invariant under the obtained realization of the 
algebra A3.5 coincides with ( |3.5|) 

whence it follows that in (4.47) e = 1 and in ( [4.48| ) 

F = —uu x + G{u x ). 
Thus the algebra ^35 = (d x , dt, td x + d u ) is the invariance algebra of the nonlinear PDE 

ih = u xx - uu x + G(u x ). 

Analysis of the cases when the operators Q2,Q2 form bases of the algebras A\ x , A\ x is 
carried out in a similar way. As a result, we get three more realizations that are invariance 
algebras of PDEs of the form fll.lQ 

Ais = (d u ,d t ,td u + Xd x ), 

u t = u xx + \~ 1 x + G(u x ), A > 0; 
^3.5 = (d u ,d x ,xd u + b(t)d x ), 

u t = u xx --b(t)u 2 x + G(t), b(t)^0; 

A t.5 = (d u ,d X} xd u + \d t ), 

u t = u xx + G(u), u = t-\u x , X ^ 0; 
Al s = (d u + 2Xtd x ,d x ,xd u + 2Xt[td t + xd x -ud u \), 

fh = u xx - 1Xuu x + r 3 G(u), u = u x t 2 - — , A ^ 0. 
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Next we consider the solvable algebras. These algebras have a common feature, namely, 
they contain commutative two-dimensional subalgebras with basis operators Qi,Q2- That is 
why, analysis of these algebras is similar to that of the algebra A 35 . 

Consider, for example, the algebra A 39 . Since the admissible pairs of the operators Qi, Q 2 
are known, all what should be done is to check the commutation relations 

[Qi, Qs] = Qu [Q 2 , Qa] = qQ2, 0<\q\< 1, (4.49) 

the operator Q 3 being of the form ( |3.6| ). 

Let the operators Qi,Q2 form a basis of the algebra A\ 1 . If Qi = d t , Q2 = 9 X , then 
it follows from commutation relations ([4.49 ) that within transformations ( |4.19| ) q = ~ and 
furthermore 

Q 3 = td t + -xd x + eud u , e = 0, 1. 

After checking the condition of invariance of equation (|1 . 1|) under the obtained realization 
of the algebra A 3 5 we see that, given the relation e = 0, the invariant PDE reads as 

u t = u xx + u 2 x G(u) 

and with e = 1 the invariant PDE takes the form 

u t = u xx + G{u), u = uu 2 x . 

Provided Q\ = d x , Q2 = d t , we get from commutation relations ( f4.49| ) that q = 2. This 
contradicts to the condition < \q\ < 1. 

Let the operators Q\, Q2 form a basis of the realization A\ A . If Q\ = d t , Q2 = d u , then 

Q3 = td t + ~xd x + qud ui < \q\ < 1. 

Provided Qi = d u , Q2 = d t , we get the following form of the operator Q 3 : 

Q 3 = qtd t + -qxd x + ud u , < \q\ < 1. 

Thus we have obtained two distinct realizations of the algebra A3. 9 

Li = (d t ,d u ,tdt + -xd x + qud u ), < \q\ < 1; 

L 2 = (d u , d t , qtd t + ~qxd x + ud u ), < \q\ < 1. 



26 



These two realizations can be unified in the following way: 

Qi = d t , Q 2 = d u , Q 3 = td t + -xd x + qud u , q^0,±l. 
The corresponding invariant equation reads 

u t = u xx + x 2{q ' 1] G(uj), uj = x l ~ 2q u x . 

At last, let us consider the case when the operators Q±, Q 2 form a basis of the realization 
A\ v This case is handled in the same way as the previous one and the results are as follows. 
We get one more realization of the algebra A3. 9 whose basis is formed by the operators 

Qi = d x , Q 2 = d u + X\t\^ 1 ~ q) d x , Q 3 = 2td t + xd x + qud u , 

where q ^ 0, ±1; A G R. The corresponding invariant equation reads 

u t = u xx - ~A(1 - q)\t\-^ 1+ ^uu x + \t\^-^G(u) 

with 

u= \t\^ l ~ q \ x . 

The remaining solvable Lie algebras are handled in an analogous way. The results on 
classification of nonlinear heat conductivity equations (|1 . 1| ) admitting the three-dimensional 
Lie algebras from the second class are summarized in Table 2, where the following notations 
are used: 

^3.3 = ( d t, tdt + -xd x , -t 2 d t - txd x + x 2 d u ), 

^3.5 = (d x ,d t ,td x + d u ), 

Al 5 = (d u ,d t ,td u + Xd x ), A>0, 

Al s = (d u ,d x ,xd u + b(t)d x ), b^O, 

Aj 5 = (d u , d x , xd u + Xd t ), A 7^ 0, 

A 5 35 = (d u + 2Xtd x ,d x ,xd u + 2Xt[td t + xd x -ud u }), A/0, 

^3.6 = ($*> &t, td t + -xd x + (it + t)d u ), 

27 



A t.6 = (dx, <9 M - - In \t\d x , 2td t + xd x + ud u ), 

A l.e = ( d u, d x ,2td t + xd x + (u + x)d u ); 

A\ % = (d u , ad x , a (a)~ 1 d t + (1 + a)xd x + [(1 — a)u + x]d u ), a = a(t), d^O 
and a 2 a + 2(d) 2 = 0; 

A\ 7 = (d t , d u , td t + -xd x + ud u ); 
A l.7 = (dx, d u ,2td t + xd x + ud u ), 
A l.s = (dt, du, td t + -xd x -ud u ), 

A i.s = (dx, d u + Xtd x ,2td t + xd x -ud u ), A e R; 

1 
Ag_ g = (<9 t , 9 X , td t + -xd x ), 

A 2 39 = (d t , d x , td t + -xd x + ud u ); 

A 3 39 = (d t , du, td t + -xd x + qud u ), q ^ 0, ±1; 

A 4 39 = (d x , d u + X\t\^^d x , 2td t + xd x + qud u ), 0<|g|<l, A e R; 

A \.w = (dx, Xtd x + d u , -\(t 2 + \- 2 )d t -\txd x + (\tu-x)d u ), A^O; 

^3.11 = (dx, ad x + d u , -(d) _1 (l + a 2 )d t + {q - a)xd x + [(a + q)u - x]d u ) , 
q > 0; a = a(t), a^O and (1 + a 2 )a = 2q(a) 2 . 



Ordinary differential equations 



a 2 a + 2{af = 0, (4.50) 

(l + a 2 )a = 2q(a) 2 . (4.51) 



can be solved by quadratures. However their general solutions are defined implicitly and 
cannot be expressed via elementary functions. 
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Table 2. Equations ( |1 . 1| ) admitting three-dimensional Lie algebras from the 

second class 



Algebra 


Function F 


A 1 

A 3.3 


\u x — x~ l u x + x~ 2 G(u), u) = 2u — xu x 


A 1 


-uu x + G(u x ) 


A 2 

^3.5 


A^ + GK), A>0, G UxUx ^0 


^3.5 


-\b{t)ul + G{t), b^O 


^3.5 


G(lu), uj = t-\u x , A^O G^^O 


^3.5 


—2Xuu x + t~ 3 G(cu), to = u x t 2 - £-, A 7^ 


A 1 


2\n\u x \G(uj), uj = x~ x u x 


A 2 


\t- l uu x + \t\~iG(u x ), 


^3.6 


\t\~5G(u>), uj = t^u 2 ,, G 7^ const, y/co 


A 4 

^3.6 


—auu x + a~ 6 exp(2a~ 1 )G'(co'), uj = u x a A — |a 3 


A 1 


G(u), uj = x~ l u x , G uu j^ 


A 2 

^S.7 


\t\--2G(u x ), G UxUx ^0 



29 



Algebra 


Function F 


A 1 


x~ 4 G(oo) } u = x 3 u x , G WUJ ^ 


A 2 


—\uu x + \t \~2G(u), u = tu x , A e R, A 2 + G^ 7^ 


A 1 

^3.9 


«^G(«), G u ^ 


A 2 

^S.9 


G(w), w = u- x ul, G w ^ 


^3.9 


x^-^GH, u; = z 1 " 2 ^, G_^0 


4 4 

^3.9 


-§A(1 - q)\t\-^+i) UUx + \t\^- 2 )G(uj), u = \t\^-i) Ux , A 2 + GL/0 


4 1 
^3. 10 


-Aim* + (t 2 + A- 2 )-lG(cj), cu = Xu x (t 2 + A" 2 ) - t, A + 


4 1 

^3. 11 


— dMMa; + (1 + a 2 )~2 exp(garctana;)G(u;), cu = 1^(1 + a 2 ) — a 



The general solution of (|4.50| ) reads as 

f a exp(-2C l )dZ = At + Ai, {A,Ai}cR, A ^ 0; 
and the general solution of ( [4.51 ) is given by the formula 

/a 
exp(— 2garctan£)d£ = At + Ai, {A, Ai}cR, A ^ 0. 

One more important remark is that the obtained realizations of three-dimensional Lie 
algebras are inequivalent. This means, in particular, that the corresponding invariant equa- 
tions are inequivalent as well. 
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5 Complete group classification of equations ( |t7l| ) in- 
variant under four-dimensional Lie algebras 

In this section we carry out group classification of nonlinear heat conductivity equations 
Q1.1D admitting four- dimensional Lie algebras. To this end, we use the known classification 
of abstract four- dimensional Lie algebras pB| . Furthermore for each invariant equation we 
compute the maximal in Lie's sense symmetry algebra thus completing the classification. 

As calculations performed for constructing inequivalent realizations of four- dimensional 
Lie algebras within the class of operators (|3.6|) are essentially the same as those used when 
we study three-dimensional ones, we will concentrate on giving the final results omitting 
calculation details. As above, we should differentiate between the cases of decomposable 
and non-decomposable four- dimensional Lie algebras. 

5.1 PDEs flTTTp invariant under decomposable four-dimensional 
Lie algebras 

The class of decomposable four-dimensional Lie algebras (regarded in a sequel as the first 
class) contains twelve algebras: 4A\ = A 3A © A u A 2 , 2 © 2A 1 = A 3 , 2 © A x , 2A 2 , 2 = A 2 , 2 © 
A 2 . 2 , A 3 _i © A\ (i = 3, 4, . . . 11). We preserve the notations of the previous section. What is 
more, A 3A = (Q u Q 2 , Q 3 ) (i = 1, 2, . . . , 11), A x = (Q 4 >. 



An analysis shows that within the class of operators ( |3.6|) there are four inequivalent 
realizations of the algebra 2A 22 which are invariance algebras of PDEs of the form 
We give these realizations below together with the corresponding invariant equations. 



(5.1) 
(5.2) 



2^2.2 " 


= (-tfy ~ 2 xd *' 5i ' 9 «» eU du), 




u t = u xx -u 2 x + —u x , A e R; 

X 


2^2.2 = 


= (-2td t - xd x , d x , d u , e u d u ), 




11 

u t = u xx u 2 x + A __, A e R; 


2^2.2 = 


= (-2td t -xd x ,d x ,-ud u + \\J\t\d x ,d u ), 




u t = u xx H t= In 1^1 + —j=\ 

A^\t\ J\t\ 




e = 1 for t > and e = -1 for t < 0, A ^ 



(5.3) 

/3GR; 
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2^2.2 



{d x -ud u ,d u ,-d t ,e xt d x ), 



Ut =M a 



Xu x (x + \n\u x \), A 7^ 0. 



(5.4) 



Next, the algebra A 3 3 © Ai has one realization which is the symmetry algebra of PDE 
belonging to the class (|1.1|) 

(d u td t + l - X d x ,-ed t -t X d x + Jd u) ®(d u ). 

What is more, the corresponding invariant equation reads as 

1 



Ut = u xx + -u x — x x u x + Ax 2 , A G R. 



(5.5) 



At last, there exists a realization of the algebra A 39 © A\ such that it is admitted by an 
equation of the form fll.lp , namely, 

(dt, d x , td t + -xd x ) © (ud u ). 

The corresponding invariant equation ( |1 . 1| ) is given below 

ik = u xx + \u~ l u 2 x , A 7^ 0. (5.6) 

All other decomposable four- dimensional algebras either have no new realizations or these 
realizations are not admitted by PDEs of the form (|1.1| ). 

Next, we carry out the complete group classification of PDEs ( |5.1D -( |5!5|) . 

Equation (5*H). 

As the equation under study contains no arbitrary functions, computing its maximal 
invariance algebra is an easy task. Performing the necessary calculations in order to solve 
( |3.5| ) yields that this algebra is infinite-dimensional. The forms of its bases operators depend 
essentially on A and are given below 



1. A ^ 0,2 





= -tdt - -xd x , X 2 = d u X 3 = 


1 


= 2t 2 d t + 2txd x + 


~\x 2 + (1 + X)t 



o„ 



o„ 



Aoo = g(t,x)e u d u , g t = g xx + -g x ; 

x 
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2. A = 0. 



3. A = 2. 



X 1 

x 4 
x 5 

Ann 



-td t - -xd x , X 2 = d t , X 3 = d u , 



2t 2 d t + 2txd x + 



-x z + t 



d u , 



td x + -xd u , X 6 = d x , 
g(t,x)e u d u , g t = g X x] 



X 1 = -td t - -jxd x , X 2 = d t , X 3 = d u , 

I 
X 4 = 2t 2 d t + 2txd x + -(x 2 + 3t)d u , 

X 5 = td x + -(x+-t) d u , X 6 = d x + -d u , 
2 V x J x 

2 
Xoo = g(t,x)e u d u , g t = g xx + -g x . 

x 

Note that the operators Xi, X 2 , X 3 and X^ with g = 1 form a basis of the algebra 2A\ 2 . 
The change of variables 

t — t, x — x, u — u — \n\x\, 

reduces the third case to the second, which means that we have two inequivalent equations 

Ut = u xx — u x ; 

u t = u xx + -u x -u 2 x , A 7^ 0,2. 



These equations are reduced to linear PDEs 



v t 



v xx + -v x , A ^ 0, 2. 
x 



with the help of the change of variables 

u = — ln|f|, u — u(t,x), v — v(t,x). 



(5.7) 
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Thus nonlinearity in equation ( |5.1| ) is not essential. 
Equation (|5.2|) 



This equation is also linearized with the aid of the change of variables (|5.7| ) to become 

A 

Vt = v xx + — j=V x . 

\t\ 



Equation (J^ 

The algebra 2A 2 2 is the maximal in Lie's sense algebra admitted by this PDE. 

Equation (|5.4|) 



Again, the algebra 1A\ 2 is the maximal symmetry algebra admitted by the equation in 
question. 

Equation (|57^ 



Performing the change of variables 

u = 41n|f|, u = u(t,x), v = v(t,x), 
yields for the function v the linear PDE 

Vt = v xx — x~ 1 v x + A\x~ 2 v . 
Equation (|5^ 



Making the change of variables 

u = ln|«|, v = v(t, x), u = u(t,x), 
reduces the equation under study to the modified Burgers equation 

v t = v xx + (A + l)v 2 x . 

The latter is locally equivalent to the linear heat conductivity equation. 

Summing up, we conclude that the class of PDEs (|1.1| ) contains only two equations ( |5.3D 
and ( |5.4|) which are essentially nonlinear and invariant under four-dimensional decomposable 
Lie algebras. And what is more, the algebras 1A\ 2 and 1A\ 2 are their maximal symmetry 
algebras. 
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5.2 PDEs (|1.1|) invariant under non-decomposable four-dimensi- 
ongal Lie algebras 

The set of inequivalent abstract four-dimensional Lie algebras contains ten real non-decom- 
posable Lie algebras A Ai = (Q 1 , Q2,Q3, Qa) (i = 1, . . . , 10) fl26f . We give below non-zero 
commutation relation determining these algebras 

A-a.i : [Q2,Qa] — Qii [Q?,iQi\ — Q2] 

Ai.2 ■■ [Qi,Qa] = qQi, [Q2,Qa] = Q2, 

^4.3 : [Qi, Qa] = Qi, [Q3, Qa] — Q2] 
A 4 a : [Qi,Qa] = Qi, [Q2,Q4] = Qi + Q2, 
[Q3,Q4] = Q2 + Q 3 ; 

-44.5 : [Ql,Q4]= Ql, [Q2,Q4}= qQ2, 

[Qs,Q4]=pQ3, -l <p<q<l, pq 7^0; 

A4.6 : [Qi, Q4] = qQi, [Q2,Q4]=pQ2-Q3, 

[Q3,Q4] = Q2+pQ3, q^o, P >o- 
A4.7 ■■ [Q2,Q3] = Qi, [Qi,Q4] = 2Qi, 

[Q 2 ,Q 4 ] = Q 2 , [Qs,Qa] = Q2 + Qs; 

^4.8 : [Q2,Q 3 ] = Qi, [Qi,Q4] = {l + q)Qu 

[Q2, Q4] = Q2, [Q3, Qa] = qQ 3 , \q\ < i; 

A4.9 : [Q2,Qs] = Qu [Qi,Q4] = 2qQu 

[Q2,Q4] = qQ2 -Q3, [Q3, Qa] = Q2 + qQ3, q>0; 

A4.10 : [Qi, Q3} = Qi, [Q2, Q3] — Q2, 

[QuQ4] = -Q2, [Q2,Q4]=Qi. 



Solving the above commutation relations within the class of operators ( |3.6| ), simplifying 
the obtained expressions for Qi, ■ ■ ■ ,Qa with the help of appropriate equivalence transforma- 
tions and solving the invariance conditions ( |3.5| ) for thus obtained operators yields that there 
are eleven realizations of non-decomposable four-dimensional Lie algebras that are symmetry 
algebras of PDEs of the form ( |1 . 1| ) . Namely, 

A\ A = (d u ,d x ,d t ,td x + xd u ); 

A\ 2 = (d t ,d u ,d x ,2td t + xd x + (u + x)d u ); 
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A l.2 = (d x ,d u ,d t ,td t + -xd x + (u + t)d u ); 

A\_ 3 = (d u , d x , d t , td t + ud u ); 

A\ 5 = (d t , d x , d u , td t + -jxd x + kud u ), k^O,-, 1; 

•4i7 = (d u ,d x ,xd u --ln\t\d x ,2td t + xd x + 2ud u ); 

A\ 8 = (d x , d t , td x + d u , td t + -xd x - -ud u ) ; 

1 3 

Ai 8 = (d u , d t , td u + \d x , td t + -xd x + -ud u ), A > 0; 

^4.8 = (d u ,d x ,xd u + X\t\^ 1 -^d x ,2td t + xd x + (l + q)ud u ), \q\^l, A ^ 0; 
^Is = (&u, d t ,xd u + \d t ,2td t + xd x + 3ud u ), A ^ 0; 

Ai 9 = (d u ,d x ,xd u + ad x ,-(a)-\l + a>)dt + ( q -a)xd x + {2 q u-^]d u ), 

where q > and the function a = a(t), a ^ is a solution of ordinary differential equation 

fUD- 

Further analysis shows that PDE ( |1.1| ) admitting the algebra A\ 1 is linearizable. All 
the remaining invariant equations are essentially nonlinear and the above algebras are their 
maximal in Lie's sense symmetry algebras. 

We present all the results on classification of inequivalent essentially nonlinear PDEs 
( |1 . 1|) that are invariant with respect to four-dimensional Lie algebras (decomposable and 
non-decomposable) in Table 3, where we use the following notations: 

2^2.2 = (-2td t -xd x ,d x ,-ud u + \\J\t\d x ,d u ), \^fr, 

2A\ 2 = (d x -ud u ,d u ,jd t ,e xt d x ), A^O; 

A l.2 = (d t ,d u ,d x ,2td t + xd x + (u + x)d u ); 

A\ 2 = (d x , d u , d t , td t + -xd x + (u + t)d u ); 

A4.3 = (d u ,d x ,d t ,td x + ud u ); 

I 1 

A\ 5 = (d t , d x , d u , td t + -xd x + kud u ), k^O,-, 1; 

Al, = (d u ,d x ,xd u - 1 -ln\t\d x ,2td t + xd x + 2ud u ); 

1 1 

A^ 8 = (d x , d t , td x + d u , td t + -xd x - -ud u ) ; 
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1 3 

A\ A = (d u , d t , td u + Xd x , td t + -xd x + -ud u ), A > 0; 

Al 8 = (d u ,d x ,xd u + \\t\^ 1 -^d x ,2td t + xd x + (l + q)ud u ), \q\^l, A ^ 0; 
At 8 = (d u ,d t ,xd u + \d t ,2td t + xd x + 3ud u ), A ^ 0; 

^4.9 = (d u , d x , xd u + ad x , -(d) _1 (l + a 2 )d t + (q - a)x9j. + [2gw - -x 2 ]<9„), 

where g > and a = a(t), a ^ is a solution of ( f4.51|) ; 
AGg(l,l) = (d x ,td x + d u ,d t ,-2td t - xd x + ud u ,t 2 d t + txd x - (tu- x)d u ). 

6 Further algebraic analysis 

In this section we prove that there are no essentially nonlinear equations of the form ( |1.1[ ) 
that admit invariance algebra of the dimension higher than 4. This means that the above 
obtained group classification of invariant PDEs (|1 . 1| ) is complete. 

Our considerations are purely algebraic and are based on the Levi-Maltsev theorem claim- 
ing that any Lie algebra over the field R or C can be decomposed into a semi-direct sum of 
a maximal solvable ideal N and semi-simple subalgebra S. This means that the problem of 
classification of abstract Lie algebras reduces to classifying 

• solvable Lie algebras, 

• semi-simple Lie algebras, 

• algebras that are semi-direct sums of semi-simple and solvable Lie algebras. 

We consider the above enumerated cases separately. 

Case 1. Solvable Lie algebras. 

As far as we know, the problem of classification of abstract solvable real Lie algebras has 



been completely solved for solvable Lie algebras of the dimension up to five (see, e.g., [f26| , |27|). 
For higher dimensional solvable Lie algebras only partial results have been obtained |3C , [31) . 
The main difficulty is that a number of non-isomorphic solvable n-dimensional Lie algebras 
increases rapidly with increasing n. For example, there are 67 types of five- dimensional 



solvable Lie algebras |27| and 99 types of six-dimensional solvable Lie algebras having a 



nilpotent element |3T| . This is why, the problem of exhaustive classification of solvable Lie 
algebras of the dimension n > 5 is a 'wild problem'. However, in the case under study it 
is possible to carry out such a classification due to the fact that we are looking for rather 
specific realizations of the solvable Lie algebras. 
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Table 3. Nonlinear PDEs ( |1.1| ) admitting four-dimensional Lie algebras 



No. 


Equation 


Maximal 

invariance 

algebra 


1 


e = 1 for t > 0, e = -1 for t < 0, (3 e R, A ^ 


2^2.2 


2 


u* = u ra - Au^rr + In u x |), A 7^ 


2^2.2 


3 


Ut = u xx + Aexp(-M x ), A 7^ 


A 1 


4 


^i = u xx + 2 In lu^l 


A 2 

A4.2 


5 


Ut = u xx — u x In \u x \ + Am x , A G R 


A 1 

A4.3 


6 


2k-2 

u t = u xx + Xu?- 1 , A ^ 0, fc / 0, |, 1 


^4.5 


7 


1 2 


A4.7 


8 


\ I I 3 
W( = M„ — WU X + A M x 2 

A^O 
A = 


A4.8 

AG|(1,1) 


9 


w* = u xx + X~ l x + mJ\u x \, A > 0, m/fl 


A 2 
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No. 


Equation 


Maximal 

invariance 

algebra 


10 


A ^ 0, \q\ ^ 1, e = 1 £ > 0, e = -1 £ < 


^4.8 


11 




^4.8 


u t = u xx + m^\t - Xu x \, X-m^0 


12 


Ut = u xx — \diu 2 x + (A — «)(1 + a 2 ) -1 , A G R 


A 1 

A4.9 



Our considerations are based on the well-known fact that for any solvable Lie algebra L n 
with dim L n = n over the field of real numbers we can construct a composition series for L n 



L C Li C • • • C L n _i C L n , 



where each algebra L iy dim Li = i (1 = 0, 1, . . . , n — 1) is an ideal in the algebra L i+i . Hence 
we easily get the following assertion. Suppose that there exist realizations A 1 , A 2 , . . . ,A 
of solvable Lie algebras within a given class of Lie vector fields V of the dimension not 
greater than m and, furthermore, realizations of the dimension m + 1 do not exist. Then the 
realizations A 1 , A 2 , ... , A N exhaust a set of all possible realizations of solvable Lie algebras 
within the class V. 

According to the results of Section 5 there are twelve realizations of solvable four- 
dimensional Lie algebras within the class of operators ( |3.6|) . If we will prove that there 
are no realizations of solvable five-dimensional Lie algebras within the class ( |3.6| ) which are 
invariance algebras of PDE of the form ( |1 . 1|) , then in view of the above assertion we conclude 
that the obtained realizations of solvable Lie algebras of the dimension n < 4 exhaust the 
set of all possible realizations of solvable Lie algebras in the case under study. 

First, we investigate the case when a five- dimensional solvable Lie algebra is a direct sum 
of four- and one-dimensional solvable Lie algebras. We consider in more detail the realization 
2Al 2 , where 



ei = -2td t - xd x , e 2 = d x , e 3 = -ud u + A J\t\d x , e 4 = d u 
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Taking the basis element e§ in the general form ([O 

e 5 = 2a(t)d t + (a(t)x + b{t))d x + f(t, x, u)d u 

and checking the commutation relations [e^es] = 0, (i = 1,2,3,4) yield that e$ = J\t\d x . 

Inserting this expression into the invariance criterion (|3.5| ) we arrive at the contradictory 

equality 

1 

-.u x = 0. 
t\ 



Consequently, the algebra A\ 2 cannot be extended to a realization of five- dimensional solv- 
able Lie algebra admitted by PDE of the form fll.ll) . The same assertion holds for the 
realizations aj\ 2 o> -^a. !•> ^*~a 8' 4 8' 4 9* 

Furthermore, the realizations A\ 2 , A% 2 , A\ 3 , A\ 5 , A\ s , A\ 8 cannot be extended to real- 



izations of five-dimensional solvable Lie algebras within the class of operators (J3T 

Next, we turn to the case of indecomposable five- dimensional solvable Lie algebras. 
According to the classification given in ||27j] there are five types of indecomposable five- 
dimensional solvable Lie algebras 

1) nilpotent algebras, 

2) algebras having one non-nilpotent basis element and containing the commuting ideal 
AA U 

3) algebras having one non-nilpotent basis element and containing the ideal A 35 © Ax, 

4) algebras having one non-nilpotent basis element and containing the ideal Ao, 

5) algebras having two nil-independent basis elements (two basis elements are called nil- 
independent if there is no linear combination of these which is nilpotent). 

Five-dimensional solvable algebras of the first type contain either a four- dimensional 
commuting radical or a radical that is isomorphic to the decomposable algebra A\ © A3.5. 
Consequently, realizations of these algebras which could be invariance algebras of PDE of 
the form ( |1 . 1|) do not exist. Similar reasonings yield the same statement for the algebras of 
the second, third and fourth types. 

Consider the algebras of the fifth type. Let ei, e 2 , e 3 , e 4 , e 5 form a basis of an algebra 
of this type. Then inequivalent abstract five-dimensional solvable Lie algebras having two 
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nil-independent basis elements read 

Li : [ei,e 4 ] = ei, [e 3 , e 4 ] = /?e 3 , [e 2 ,e 5 ] = e 2 , 

[e 3 ,e 5 ]=7e 3 , /3 2 + 7 V 0; 
L 2 : [ei,e 4 ]=aei, [e 2 ,e 4 ] = e 2 , [e 3 ,e 4 ]=e 3 , 

[ei,e 5 ] = ei, [e 3 ,e 5 ] = e 2 ; 
L3 ■ [ei,e 4 \=ae 1 , [e 2 ,e 4 ] = e 2 , [e 3 ,e 4 ]=e 3 , 

[ei,e 5 ] = 6ei, [e 2 ,e 5 ] = -e 3 , [e 3 ,e 5 ] = e 2 , a 2 + 5 2 ^ 0; 

£4 : [e 2 ,e 3 ]=e x , [ei,e 4 ] = ei, [e 2 ,e 4 ] = e 2 , 

[e 2 ,e 5 ] = -e 2 , [e 3 ,e 5 ] = e 3 ; 

£5 : [e 2 ,e 3 ]=ei, [ei,e 4 ] = 2ei, [e 2 ,e 4 ]=e 2 , 

[e 3 , e 4 ] = e 3 , [e 2 , e 5 ] = -e 3 , [e 3 , e 5 ] = e 2 ; 

^6 : [ei,e 4 ]=ei, [e 2 ,e 5 ] = e 2 , [e 4 ,e 5 ] = e 3 ; 

L 7 : [ei,e 4 ] = e t , [e 2 ,e 4 ] = e 2 , [ei,e 5 ] = -e 2 , 

[e 2 ,e 5 ] = ei, [e 4 , e 5 ] = e 3 . 

Note that we give non-zero commutation relations only. 

The algebra L\ contains a radical isomorphic to the decomposable four-dimensional Lie 
algebra A\®Az$. Next, the algebra Lq contains a radical isomorphic to the algebra 2/li©/l 2 . 2 . 
At last, the algebra L 7 contains a radical isomorphic to the algebra A\ © A3.7. Hence we 
conclude that there are no realizations of the algebras L\,L§, L7 which are invariance algebras 



of PDE (1TT 



The algebra A\ 5 gives a realization of the algebra A4.5 with q = |,p 7^ 0, |, 1. On the 
other hand, the algebras L 2 , L 3 contain a radical isomorphic to the algebra A 4 .5 with g = 1. 
Hence it follows that the realization A\ 5 cannot be extended to yield a realization of the 
algebras L 2 ,L 3 . 

The algebra L 4 contains a radical isomorphic to the algebra A 48 with q = 0. To meet 
this requirement we have to choose for this radical the realization A\% with q = 0, namely, 

ei = d u , e 2 = <9 X , e 3 = x^ + A|t| *d x , e 4 = 2t<9 t + xd x + ud u . 



Checking commutation relations for an operator e§ of the form (|3.6|) shows that the realiza- 
tion y4| 5 cannot be extended to give a realization of the five-dimensional algebra L 4 . 

The algebra L5 contains a radical isomorphic to the algebra A^s (q = 1). However, there 
are no realizations of the algebra y4 4 .g which might yield a realization of this radical. Hence 
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we conclude that there are no realizations of the algebra L5 within the class of operators 



(PD- 

Case 2. Semi-simple Lie algebras. 

As proved by Cartan, any real or complex semi-simple Lie algebra is decomposed into 
a direct sum of mutually orthogonal simple algebras. In view of this fact, the problem of 
classification of abstract semi-simple Lie algebras reduces to classifying simple Lie algebras 



see, e.g. [J32|). The classification of simple Lie algebras is well-known. There are four series 



of non-exceptional complex simple Lie algebras A n , B n , C n , D n and five types of exceptional 
Lie algebras. 

The lower dimensional semi-simple Lie algebras are connected by the following isomor- 
phisms [ |32]j : 

su{2) ~ so(3) ~ sp{l), sl(2, R) ~ su(l, 1) ~ so(2, 1) ~ sp(l, R), 

so(5)~sp(2), so(3,2) ~ sp(2,R), so(4, 1) ~ sp(l, 2), 

so(4) ~so(3)©so(3), so(2,2) ~ sZ(2,R) © sZ(2,R), (6.1) 

sZ(2,C) ~so(3, 1), sw(4)~so(6), s/(4, R) ~ so(3, 3), 

sw(2,2) ~so(4,2), su(3, 1) ~ so*(6), sm*(4) ~ so(5, 1). 

It turns out that si (2, R) ~ sm(1, 1) ~ so(2, 1) ~ sp(l, R) are the only real forms of the 
algebras given in ( |6.1| ) that have realizations within the class of operators (|3.6|) . The reason is 
that all other algebras contain the subalgebra so(3) and the latter has no realizations within 



the class ( |3.6j ). Next, all the real forms of higher dimensional non-exceptional simple Lie 
algebras contain the algebra so(3) as a subalgebra. Consequently, they have no realizations 
within the class of operators (|3.6|). 

The exceptional simple Lie algebras have no realizations within the class of differential 
operators of the form ( |3.6| ) . 

Consequently, the only semi-simple algebras that might be admitted by PDE of the form 
( |1 . 1|) are algebras of the form 

s/(2,R), s/(2,R)©sZ(2,R), sl{2, R) © sl(2, R) © sl{2, R), . . . 

As straightforward calculation shows, there are no PDEs of the form (|1 . 1|) invariant with 
respect to the algebra sl(2, R) ©s/(2, R). Hence it follows, that the only semi-simple algebra 
that might be admitted by (|1.1|) is the three-dimensional algebra s/(2,R). 

Case 3. Semi-direct sums of semi-simple and solvable algebras. 
The algebras of the type considered can be split into two classes. 
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• algebras which are decomposable into direct sums of semi-simple and solvable algebras, 

• algebras which cannot be decomposed into direct sums of semi-simple and solvable 
algebras. 

As shown above, there exists only one realization A\ 3 of a semi-simple algebra which is 



an invariance algebra of an equation of the form ( |1 . 1| ) . It is a realization of simple algebra 
A3. 3 isomorphic to the algebra sZ(2, R). If we will try to extend this realization to get a 
realization of a direct sum of semi-simple and solvable Lie algebras, then we will have to 
stop at the first step, since the realization A\ 3 © A\ is an invariance algebra of linear PDE 
(see Section 5). 

Turn now to the algebras which are not decomposable into a direct sum of semi-simple 
and solvable Lie algebras. According to the above results of the previous two cases, their 



dimension cannot be higher than 3 + 4 = 7. In the paper |S3| a complete classification of the 



algebras which are semi-direct sums of semi-simple and solvable Lie algebras and have the 
dimension n < 8 is obtained. Analysis of these algebras shows that they have no realizations 
within the class of operators (|3.6|) that are invariance algebras of PDE of the form (|1 . 1|) . 
Summing up we conclude that there are no real Lie algebras of the dimension n > 5 



which are invariance algebras of essentially nonlinear PDEs belonging to the class (|Ll 
This means that our classification of nonlinear PDEs ( |1 . 1| ) invariant under the one-, two-, 
three- and four- dimensional Lie algebras gives the complete description of heat equations 
( |1.1|) possessing non-trivial Lie symmetries. 

7 Comparison to other classifications 

Here we briefly review the earlier results on classification of invariant PDEs belonging to the 
class ( |1 . 1| ) . We will show that all of them can be derived from equations given in Tables 1-3 
(either directly or via local transformations of dependent and independent variables). 

The problem of group classification of the nonlinear heat conductivity equation with a 
nonlinear convection term 

u t = [K(u)u x ] x + [®{u)] x (7.1) 

has been considered in [22 , ^3]] . Evidently, provided K(u) = 1, it is included into the class 



(ED- 

Next, Dorodnitsyn |Hj] has classified invariant nonlinear heat conductivity equations with 



nonlinear source 



dT d ( dT\ 
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Again, this equation with K(u) = 1 belongs to the class ( |1 . 1| ) . Note that an analogous 
problem for the two- and three-dimensional PDEs of the type ( |7.2|) has been solved in [29 . 
The papers [ES] are devoted to symmetry analysis of nonlinear PDEs of the form 

u t = [A(u)u x } x + B(u)u x + C(u). (7.3) 



Nonlinear PDE ( \l.3\) is a natural generalization of equations (|7.1| ), ( |7.2|) and, furthermore 



is contained in the class of PDEs ( |1.1| ) provided A(u) = 1. 



Gandarias [24| has carried out group classification of equation 

u t = {u n ) xx + g{x)u m + f(x)u"u x , n^O (7.4) 



that is also included into the class ( |1 . 1|) , provided the condition n = 1 holds. 

7.1 Group analysis of equation ( |7 . 1| ) 

According to [^2|, ^3] the results on group classification of equation ( |7.1| ) under K{u) = 1, 
namely of equation 

u t = u xx + [$(w)U (7.5) 

can be summarized as follows. The maximal invariance algebra admitted by PDE ( [7.5| ) 
under an arbitrary function $, ^ ^ is the two-dimensional Lie algebra (d t , d x ). Extension 
of the invariance algebra is only possible, provided 

1) $> = (3u u , Q ncw = 2{1 - u)td t + {1 - u)xd x + ud u , 

2) $ = (3 In u, Q new = 2td t + xd x + ud u ; 

3) $ = /3e^, Q^ = td t + \xd x -^-d u . 

Here v ^ 0, 1, 2, /3 e R. 

Note that equation (|7.5| ) with $ = /3w^, where z/ = 2 coincides with the Burgers equation 
which maximal symmetry algebra is five-dimensional and is isomorphic to the full Galilei 
algebra. 

For the first case ($ = j3u u ) the invariance algebra is isomorphic to the algebra A3. 9 
(q = |). This isomorphism is established by choosing the basis operators as follows 

Qi = d u Q 2 = d x , Q 3 = td t + -xd x + ud u . 
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Furthermore the change of variables 

t = t, x = x, v = u 2{l ~ v) 

transforms the above realization to become A 2 9 . With this transformation the corresponding 
invariant equation ( [T.5|) takes the form 

2v — 1 _-, 2 2v _i 

«t = " ra + 7771 V w ^ + P^T\ \ v * u *> 7 - 6 

2(1 — v) 2(1 — v) 

which is a particular case of the equation invariant with respect to the algebra A\ 9 from 
Table 2. 

Given the condition $ = (3lnu, the invariance algebra of ( |7.5|) is also isomorphic to the 
algebra A3.9 (g = |), its basis being chosen in the following way: 

1 1 

d t , d x , td t + -xd x + -ud u . 



The change of variables 



t = t, x = x, v = u 2 



reduces the corresponding invariant equation ( |7.5|) to the form 

v t = v xx - ^ l v 2 x + pv~2 Vx . 

The latter is, evidently, a particular case of PDE invariant with respect to the algebra A\ 9 
from Table 2. 

At last, for the third case the invariance algebra is also isomorphic to the algebra A 39 (q = 
|) and is reduced to the realization A\ 9 with the help of the change of variables t = t, x = 
x, v = e~ 2uu . The corresponding invariant equation (|7. 5| ) with this change of variables takes 
the form 

v t = v xx - v~ x vl + (3vv~iv x , 

which is a particular case of PDE invariant with respect to the algebra A\ 9 from Table 2. 
Summing up we conclude that the group classification of PDE (7.5) within the equivalence 



relation follows from our classification of equations invariant under the Lie algebra A\ 9 if 
we put in these 

G{uS) = XiU + A 2 a/cJ, uj = u~ 1 u 2 ,, {Ai,A 2 }cR. 
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7.2 Group analysis of equation ( |7.2| ) 

The results on group classification of ( |7.2| ) with K{T) = 1, namely for PDE of the form 

u t = u xx + F(u) (7.7) 

given in [^TJ can be formulated in the following way. Provided the function F, ^Jr 7^ is 
arbitrary, the maximal invariance algebra of ( [7.71) is the two-dimensional Lie algebra (dt, d x ). 
Extension of the invariance algebra is only possible provided 

1) F = ±e\ Q ncw = td t + ^xd x -d u , 

2) F = ±u n , Q new = td t + -xd x -ud u ; 

2 n — 1 

3) F = 5ulnu, 6 = ±1, Ql w = e 5t [d x - 6 -xud u ], Q 2 new = e 5t ud u . 

Note that the classification results yielding linear invariant PDEs are neglected here. 
Consider first the case 1 . Then the change of variables 

t — t, x = x, v = e~ u 

reduces the invariance algebra to become A\ 5 and, furthermore, the corresponding invariant 
equation ( |7.7|) takes the form 

v t = v xx -v^vIt 1. 

For the second case, there is the change of variables 

t = t, x = x, v = u l ~ n , n^l 

that reduces the invariance algebra to become A| 5 . The corresponding invariant equation 

(|7.7|) takes the form 

n -i 2 ' 

v t = v xx H -v v x ± 



n - 1 x 1 - n 

The above two PDEs are particular cases of the equation invariant with respect to the 
algebra A| 5 from Table 2. 

At last, in the third case the maximal invariance algebra is four- dimensional. Making 
use of the change of variables 

T = - 6 -e~ 2St , Z = e~ 5t x, v = e- St [\n\u\ + S -x 2 } 
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we become convinced of the fact that the invariance algebra is equivalent to A\ 8 with q 
0, A = 2S\/2. The corresponding invariant equation (|7.7|) is reduced to the form 



eS , ,_i / , 2 e5 



U «--H T I *W) + 



2 i i v u 2 ^i i , 
where e — 1 for r > and e — —1 for r < 0. Making the second change of variables 

yields the equation under the number 9 from Table 3 with A = 2y / 2<5, q = and 



Similar analysis of classification results for PDEs ( [7.3|) | 2q| and ( |7.4|) [p4|1 shows that all 
the invariant equations obtained there can be derived from invariant PDEs given in Tables 
2, 3 under appropriate changes of variables. We unable to present here the corresponding 
calculations in a compact form, since they are extremely lengthy (just a precise formulation 
of classification results obtained in |25|, [24]] requires several pages, to say nothing of a space 
needed to give a detailed analysis of these). 

8 Concluding Remarks 

We have carried out group classification of nonlinear heat transfer equations of the form 
(|1.1|) and proved that essentially nonlinear PDEs ( |1 . 1| ) admit at most four-parameter in- 
variance group. Furthermore, we have established that there are three classes of equations 
( |1 . 1| ) invariant with respect to one-parameter groups (formulae ( |4.4| )-( |4~oT) ), seven classes 
of equations (|1.1|) invariant with respect to two-parameter groups (formulae ( |4.10 )-( 4.12| ), 



Q4.14| ), ( |4.15| ), ( |4.17| ), ( [4.1^ )), twenty eight classes of equations ( [Lip invariant with respect 



to three-parameter groups (Tables 1, 2) and twelve classes of equations ( |1 . 1| ) invariant with 
respect to four-parameter groups (Table 3). 

We concentrate on studying essentially nonlinear heat conductivity equations since the 
linear case is well investigated. However, it is fairly simple to recover the corresponding 
results within the framework of our approach. Consider the most general linear PDE of the 
parabolic type in one spatial variable 

u t = f(t, x)u xx + g(t, x)u x + h(t, x)u. (8.1) 
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The most general infinitesimal operator of the symmetry group admitted by ( |3.1| ) reads 



as 



Q = T(t)d t + X(t, x)d x + (U(t, x)u + u (t, x))d u , 
where T,X,U are arbitrary smooth functions and Uq is an arbitrary solution of (|8.1|). As 



usual, we neglect the trivial symmetry u (t,x)d u and put u = 0. Next, the equivalence 
group of the class of PDEs ( |S.lD has the form 

t = F(t), x = G(t, x), u = H(t, x)u. 

Using these facts it is straightforward to check that the list of inequivalent one-dimensional 
Lie algebras admitted by ( ft.l|) is exhausted by the following three algebras 

Ai = (0 s ), A 2 = (dt), A 3 = (U(t,x)d u ). 

As equation (|8.1|) is linear, it admits the one-dimensional Lie algebra ud u with arbitrary 
f,g,h. Consequently, any two-dimensional algebra is reduced to the one of three possible 
inequivalent forms (ud u ) © A4 (i — 1, 2, 3). 

If equation (|8.1|) is invariant with respect to the algebra (ud u ) © A\, then its coefficients 
are independent of x. Hence we easily get that it is reduced to the standard heat transfer 
equation 

u t = u xx . (8.2) 



Turn next to the case of the algebra (ud u ) © A 2 . Now the coefficients of (|8.1|) are inde- 
pendent of t and, therefore, this equation can be reduced to become 

Ut = u xx + V(x)u (8.3) 

with an arbitrary smooth function V . It is a common knowledge that the above PDE has a 
symmetry algebra of the dimension higher than 2 if and only if 

V(x) = ^ + X x x 2 + X 2 x + A 3 , (8.4) 

x z 

where Ao, . . . , A3 are arbitrary constants with A0A2 = 0. Furthermore, provided Ao = 0, PDE 
(|yD, (0) is equivalent to the heat transfer equation (fT2|). If, A ^ 0, then PDE flOD , 
reduces to the following equation: 



Ut = u xx + —u (8.5) 



x 2 
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which is invariant under the four-dimensional Lie algebra 

( d t , 2td t + xd x , t 2 d t + txd x - - + - ud u , ud u ) . 

Summing up we conclude that there are three inequivalent classes of PDEs ( |S.1| ) whose 
symmetry algebras have the dimensions higher than one, namely, the heat transfer equation 
( |8.2| ) admitting the six- dimensional Lie algebra, equation ( |8.5| ) invariant with respect to 
the four- dimensional algebra and equation (|8.3|) that admits the two-dimensional algebra 
(d t ,ud u ). This completes group classification of heat transfer equations ( |1.1[ ) admitting 
nontrivial Lie symmetry. 

When classifying invariant equations ( |1.1|) we utilize as equivalence transformations local 
transformations of dependent and independent variables. Using non-local transformations, 
on the one hand, may result in reduction of equivalence classes and, on the other hand, may 
yield so-called quasi-local symmetries (for more detail on on quasi-local symmetries see, e.g. 



[01 )• Consider, as an example, the following subclass of PDEs of the form ( |1.1|) : 



u t = u xx + fi(t)u + f 2 {t,x,u x ) (8.6) 

with arbitrary smooth functions fi, f 2 - If we differentiate ( |1.1|) with respect to x and make 
a change of the dependent variable 

u x (t,x) — > v(t,x), (8.7) 

then we get a subclass of quasi-linear PDEs of the form (|1 . 1|) 

v t = v xx + fi(t)v + f 2x (t,x,v) + f 2v (t,x,v)v x . (8.8) 



Evidently, the above two classes of PDEs (|8.6|) and (|8.7|) are inequivalent in the sense of the 
definition given in Section 3, since transformation ( fj.7| ) is not local. 

The technique developed in the present paper can be efficiently applied to carry out 
group classification of arbitrary classes of PDEs in two independent variables, since their 
maximal symmetry algebras are, as a rule, low dimensional and we can use the classification 
of abstract low dimensional Lie algebras. 

These and the related problems are under study now and the results will be reported in 
our future publications. 
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